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CONVERSION  FACTORS,  NON-SI  TO  SI  UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as  follows: 


Multiply 

By 

To  Obtain 

inches 

2.54 

millimetres 

kips  (force)  per  square  inch 

6.894757 

megapascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

I.  INTRODUCTION 


1.1  Review  of  Literature  on  Wavefront  Propagation 

One  possible  way  of  studying  wave  propagation  is  to  follow  the  evolution  of 
wavefronts  in  spacetime.  A  need  to  approach  many  problems  in  that  way  arises  in  cases 
of  intensive  dynamic  loads  acting  on  material  bodies.  The  researches  are  typically  based 
on  either  the  theory  of  characteristics  or  the  theory  of  propagating  singular  surfaces.  The 
first  of  these  relies  on  the  methods  of  solution  of  hyperbolic  partial  differential  equations 
and  borrows  heavily  from  developments  in  gas  dynamics.  This  approach,  in  which  the 
first  papers  appeared  in  the  forties,  developed  rapidly  in  the  sixties,  and  is  described  in 
the  book  of  Cristescu  (1967).  This  book  is  devoted  to  the  transient  response  of 
continuous  plastic  bodies  and  remains  a  classic  reference  to  this  day.  Many  advances 
were  made  later,  particularly  in  the  class  of  one-dimensional  (1-D)  problems,  but  also  in 
various  2-D  and  3-D  axisymmetric  problems.  Various  constitutive  laws  were  considered 
-  such  as  elastic-plastic  and  elastic/viscoplastic  -  since  the  effort  was  on  both  solids  and 
soils.  This  progress  is  evidenced  in  the  monograph  by  Nowacki  (1978).  Another 
monograph  of  similar  quality  is  the  one  due  to  Wlodarczyk  (1986).  We  conclude  this 
very  brief  account  of  this  class  of  studies  on  wave  propagation  in  elastic-inelastic 
materials  with  an  observation  that  the  material  spatial  randomness  -  and  hence, 
randomness  in  constitutive  response  and  its  effect  on  wave  passage  -  has  not  been 
accounted  for. 

The  theory  of  propagating  singular  surfaces  dates  back  to  Hadamard  (1903)  who 
provided  the  celebrated  Hadamard  compatibility  condition.  This  is  a  condition  that  all 
surfaces  of  discontinuity,  strong  or  weak,  have  to  satisfy.  Evolution  of  surfaces  of 
discontinuity  within  the  framework  of  modem  continuum  mechanics  was  studied 
rigorously  by  Truesdell  and  Toupin  (1960)  and  Truesdell  and  Noll  (1965).  In  the 
following  years  this  approach  was  utilized  to  examine  the  growth  and  decay  behavior  of 
shock  waves  and  acceleration  waves  in  various  classes  of  materials.  Many  results  in  this 
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area  were  treated  in  the  book  by  Chen  (1976).  Again,  a  reader  interested  in  the  effect  of 
material  randomness  on  the  propagation  of  these  discontinuity  waves  will  come  away 
disappointed:  the  deterministic  continuum  theory  is  well  advanced,  but  the  stochastic 
problems  are  not  yet  treated.  This  statement  is  meant  not  as  a  criticism  but  only  as  an 
observation;  the  deterministic  theories  always  provide  a  basis  for  extensions  to  stochastic 
phenomena  in  physics  and  mechanics.  Since  such  extensions  are  usually  not 
straightforward,  various  existing  deterministic  approaches  have  to  be  learned  in  order  to 
find  the  best  route.  Thus,  in  the  remainder  of  this  section  we  give  short  accounts  of 
several  other  wavefront  propagation  analyses  we  considered. 

Of  special  interest  to  us  are  theoretical  researches  on  transient  waves  in 
inhomogeneous  media.  One  of  the  first  papers  in  this  category  was  due  to  Sternberg  and 
Chakravorty  (1959),  who  investigated  the  propagation  of  rotary  shear  waves  in  an 
inhomogeneous  isotropic  elastic  plate  from  a  circular  opening.  The  inhomogeneity  was 
considered  in  the  form  of  a  power-like  dependence  of  the  shear  modulus  on  the  radial 
distance  from  the  center  of  the  hole.  Indeed,  this  problem  provided  the  basis  for  other 
theoretical  techniques  and  more  general  assumptions.  Thus,  while  the  original  problem 
could  be  tackled  by  a  Laplace  transform,  a  wider  class  of  inhomogeneities  was  treated  by 
a  numerical  scheme  based  on  the  theory  of  characteristics  (Chou  and  Schaller,  1966). 
Following  this,  a  new  approach  was  developed  by  Achenbach  and  Reddy  (1967)  which 
relied  on  the  Taylor  series  representation  of  the  field  quantities  in  the  vicinity  of  the 
wavefront.  This  method  was  utilized  by  several  authors,  and  was  critically  assessed  by 
McCarthy  (1975).  In  view  of  the  drawback  of  the  method  in  the  setting  of  a 
deterministic  medium,  this  line  of  approach  does  not  seem  promising  for  media  with 
random  microstructures. 

Recently,  a  new  method  of  solving  linear  wave  equations  has  been  developed  by 
Seymour  and  Varley  (1989).  While  the  method  can  be  applied  to  problems  involving 
very  general  spatial  inhomogeneities  of  the  medium,  its  usefulness  for  us  is  restricted  by 
the  assumption  of  a  linear  elastic  constitutive  behavior. 
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Strong  material  spatial  randomness  of  geologic  media  was  the  motivation  for 
Rohani  (1982)  who  initiated  a  program  of  research  aimed  at  bringing  out  the 
microstructural  nondeterministic  effects.  In  that  study  the  problem  was  considered  as 
one  of  wave  propagation  in  a  homogeneous  medium  whose  defining  physical  constants 
are  random  variables.  In  addition,  the  characteristics  of  airblast  were  considered  random 
too.  Several  sensitivity  analyses  were  conducted  to  assess  the  relative  effects  of  these 
random  quantities  on  the  scatter  of  output  quantities.  Further  studies  in  this  vein  were 
reported  in  (Rohani  and  Cargile,  1984). 

An  attempt  to  introduce  the  spatial  variability  of  the  medium’s  properties  was 
undertaken  in  (Sadd,  Hossain  and  Rohani,  1986).  This  study  was  based  on  a  distributed 
body  concept  advanced  by  Goodman  and  Cowin  (1972)  and  the  associated  wave 
propagation  studies  conducted  by  Nunziato  and  Walsh  (1978)  and  Nunziato,  Kennedy 
and  Walsh  (1978).  The  wave  propagation  calculations  were  based  on  the  Bernoulli 
equation  for  the  evolution  of  an  acceleration  wave,  in  which  the  two  material  coefficients 
were  related  to  the  microstructure  in  the  distributed  body  concept  Similar  as  before, 
these  coefficients  were  taken  as  space  homogeneous  constants  assigned  random  values. 

1.2  Motivation  and  Scope  of  the  Present  Research 

Randomness  of  physical  and  geometrical  microscale  properties  is  an  inherent 
characteristic  of  most  materials.  It  is  definitely  the  case  with  granular  media,  and  this 
aspect  has  to  be  taken  into  account  in  any  model  aimed  at  a  more  precise  description  of 
mechanical  phenomena.  The  phenomenon  under  study  in  this  report  is  one  of  transient 
wave  propagation  in  a  granular  microstructure;  see  Fig.  la).  We  follow  the  theoretical 
approach  of  our  earlier  works  (Ostoja-Starzewski,  1984  and  1989a).  The  model  is 
developed  in  detail  in  a  one  dimensional  (1-D)  setting,  that  is,  for  a  microstructure 
represented  by  a  sequence  of  grains  whose  properties  -  mass  density  and  constitutive 
moduli  -  are  randomly  assigned  but  constant  throughout  the  domain  of  the  grain.  The 
approach  adopted  here  is  based  on  the  recognition  of  the  Markovian  property  of 
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disturbance  propagation  in  so  far  as  the  ray  kinematics  as  well  as  the  amplitude 
modulation  of  the  wave  transmitted  forward  are  concerned.  Descriptions  of  both  these 
processes  in  terms  of  Fokker-Planck  equations  are  obtained,  and  these  lead  in  turn  to  the 
formulas  that  describe  the  dependencies  of  first  and  second  moments  on  the  medium’s 
statistics  and  the  propagation  times.  The  ensuing  model  is  thus  based  on  a  semi-group 
property  of  a  Markovian  propagator,  which  reflects  a  stochastic  form  of  Huygens’  minor 
principle.  Since  this  model  is  applicable  to  any  characteristics  of  the  forward 
propagating  (transmitted)  waves  up  to  the  point  of  intersection  with  another 
characteristic,  it  can  be  used  to  investigate  1-D  problems  of  transient  response  of 
microstructures  with  various  constitutive  laws.  Two  cases  are  treated  here: 
microstructure  with  linear-elastic  grains,  and  microstructure  with  linear-hysteretic  grains. 

In  the  last  chapter  of  this  report  we  show  how  our  theoretical  approach  can  be 
extended  to  two  and  three  dimensions.  Starting  from  the  Markov  property  of  disturbance 
propagation  we  present  a  generalization  of  the  diffusion  approximation,  which,  in 
conjunction  with  solutions  of  transient  wave  problems  in  deterministic  media,  can  be 
used  to  find  the  responses  of  randomly  heterogeneous  media.  Finally,  we  introduce  the 
concept  of  local  averaging  of  wavefronts. 
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II.  1-D  WAVEFRONT  PROPAGATION 


2.1  Basic  Model 


In  this  section  we  consider  wavefront  propagation  in  a  1-D  model  of  a  granular 
medium.  That  is,  by  a  granular  medium  we  understand  a  family  of  semi-infinite  (Xi  ^  0) 
sequences  of  grains;  response  of  each  sequence  is  independent  of  that  of  others.  Each 
grain  is  taken  as  a  homogeneous  continuum  of  physical  properties,  such  as  mass  density 
p(oi)  and  elastic  modulus  E(co),  sampled  from  certain  probability  distributions.  In 
addition,  the  lengths  l(co)  of  grains  may  also  be  random.  It  is  seen  that  a  single  sequence 
of  grains  is  a  deterministic  medium  B(co),  while  a  family  B  =  (B(co);  co  e  ft)  is  a  random 
medium,  whereby  ft  denotes  an  underlying  sample  space;  see  also  Sobczyk  (1986). 

We  consider  a  space-time  graph  of  a  disturbance  propagating  in  a  class  of  media  B 
in  which  grains  are  linear  elastic  and  grain  boundaries  are  possibly  dissipative.  By  a 
disturbance  we  understand  any  single  point  of  the  pulse  p(t)  shown  on  the  t-axis  of  Fig.  2; 
this  pulse  results  in  a  wavefront  moving  in  the  material  domain. 

We  observe  that  due  to  random  speeds  of  propagation,  any  initial  disturbance,  from 
to  an  arbitrary  point  p(to)  of  the  pulse,  diffuses  in  the  space-time  within  a  forward 
causality  cone  C+.  The  cone  C+  is  a  set  of  all  possible  stochastic  ray  paths 
(characteristics)  that  originated  from  X=Xj  =0att  =  to,  which  in  view  of  the  model’s 
microstructure  are  continuous  piecewise  linear  lines.  In  the  following  we  take,  without 
loss  of  generality,  t„  =  0.  Locally,  the  gradient  of  the  characteristic  is  the  phase  velocity 


c(co)  = 


da(e,  to) 
d£ 


/p(co) 


2 


(1.1) 


Now,  we  introduce  a  dispersion  distance  defined  by 

S(t,  co)  =  X  (t,  co)  -  <X(t)>  ( 1 .2) 
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where  <X(t)>  =  cot  is  the  position  of  the  average*  path  defined  as  a  path  in  an  average 
medium.  It  follows  that  £(t)  is  a  random  walk  about  the  origin  £  =  0. 

A  model  fully  equivalent  to  this  one  may  be  formulated  by  introducing  a 
dispersion  time: 

x(X,  fo)  =  t(X,  co)-T(X)  (1.3) 

which  describes  the  Markovian  spreading  of  paths  about  the  average  path,  see  Ostoja- 
Starzewski  (1984).  In  the  latter  reference  x  was  parametrized  by  the  average  time 
T=  X/<c>,  where  <c>  is  the  ensemble  average  of  phase  velocities  c  of  the  grains. 

Let  us  define: 

'Fq0  (t)  —  set  of  all  materials  points  in  B  disturbed  in  time  £  t; 

(t)  =  3Yqo  (t)  —  set  of  all  material  points  in  B  reached  at  time  L 
Thus  Oqo  (t)  is  a  disturbance  at  t  due  to  a  cause  at  point  q0  =  (X  =  0,  to  =  0);  that  is, 
<Dqo(t)  is  a  zone  of  finite  rather  than  zero  thickness  as  would  be  the  case  of  a  disturbance 
propagation  in  a  deterministic  medium.  We  note  that  ^(t)  and  C^ft)  may  also  be 
parameterized  by  the  average  time  T  as  indicated  in  Fig.  2  -  then  we  have  ^(t)  and 

The  evolution  of  4^(0  is  determined  by  the  dynamics  of  the  £t  process  (definition 
(1.2)),  or  equivalently,  by  the  Xr  process  (definition  (1.3)).  Let  us  consider  the  process  £t 
from  now  on.  The  physical  properties  p(co)  and  constitutive  law  o(e,co)  are,  in  general, 
correlated  in  spatial  domain.  Without  loss  of  generality  we  assume  them  to  be 
describable  by  a  Markov  process  (p,a)t.  It  follows  then  from  (1.1)  that  £t  is  Markov 
only  in  the  sense  of  the  vector  process  (p,c,£)t  being  Markov.  For  simplicity,  however, 

*  The  symbol  <•>  is  used  to  denote  ensemble  averages. 
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we  shall  assume  henceforth  that  the  statistics  of  the  physical  properties  are  spatially 
homogeneous  in  the  strict  sense.  In  that  case  alone  is  Markov  and  its  transition 
function  is  time-homogeneous 

P(t,x,t  +  At,F)  =  P(At,x,F)  =  P{£(At)eF  |  ^(0)  =  x)  (1.4) 

We  thus  specify  this  process  in  terms  of  its  one-step  transition  function,  where 

At  =  <l/c>  (1.5) 


is  the  basis  of  a  (yet  to  come)  continuous  diffusion  process  approximation  of  the  discrete 
physical  process.  Also,  F  is  any  set  in  the  range  (state  space)  E  of  the  variable 
Furthermore,  since 


p(At,x,y)  =  P  (£(At) * y  1 4(0)  =  x)  =  P{£(At)  =  y-x  |5(0)  =  0)  d-6) 


the  transition  function  is  also  homogeneous  in  the  state  space.  In  the  above  we  have  used 
a  transition  density 


p(At,x,y) 


9P(At,x,y) 

dy 


(1.7) 


Now,  we  note  that  the  above  function  can  be  determined  in  terms  of  the  medium’s 
statistics 


p(At,x,y)  =  P{(o:  X(At,oo)  -  <c>  At  =  y -  x},  X(At,o>)  =  c(co)At  (1.8) 

The  dynamics  of  the  process  is  described  by  a  Chapman-Kolmogorov  equation 

p(ti  + t2,  xo,  x2)  =  Jp(ti,  xo,  Xj)p(t2,  xlf  x2)dx!  ,  (1.9) 

where  the  integration  is  over  C+(ti)  Pi  C"(t2),  and  C_(t2)  is  a  backward  dependence 
cone  for  a  time  difference  t2. 
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It  follows  from  the  theory  of  Markov  processes  that,  by  introducing  an  operator 


U(t)[f(x)]  =  j  f(y)P(t,  x,  dy),  f  e  Li  (E)  (UO) 

relation  (1.9)  leads  to  a  semi-group  property: 

U(t1+t2)  =  U(t1)U(t2)  (1.11) 

which  represents  a  stochastic  form  of  Huygens’  minor  principle  for  the  arrival  of 
disturbance  0^(0. 

At  this  stage  we  introduce  a  diffusion  approximation  of  the  £t  process,  which  is 
expressed  by  the  following: 

i)  the  Fokker-Planck  equation* 

•J^-p(t,  x)  =  -  -|-[A5(x)p(t,  x)]  +  y-^-[B^(x)p(t,  x)]  (1-12) 


in  which 

ii)  p(t,  x)  is  the  probability  density 

p(t,  x)  =  |p(t,  x',  x)dx'  (1.13) 

<5 

iii)  satisfying  the  initial  condition 

p(0,  x)  =  6(x  -  Xo)  (1.14) 

iv)  and  p(t,  x  ,  x)  being  such  that 


* 


The  time  derivative  here  is  meant  in  the  sense  of  forward  time  differencing,  i.e.  -f-f(t)  =  — +Atj — — 

a  t  At 

where  f(t)  is  an  arbitrary  function. 
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A^x) = -4— <^<<a*>)  -  m  i  m = x>  a-is) 

<ax> 

B^(x)  =  -i—  <£(<ax>)  -  £(0)]2 1  5(0)  =  x>  (1.16) 

<ax> 

-4-<{[4«°x> -  W  14(0)  =  X> S 0,  n  >  2  (1.17) 

<°x> 

Condition  (1.17)  represents  a  constraint  on  strength  of  randomness  of  the  medium’s 
microstructure  for  the  validity  of  the  diffusion  approximation. 

The  wave  amplitude  £  evolves  as  a  Markov  process  with  discontinuous  sample 
paths.  In  view  of  the  spatial  homogeneity  of  the  medium’s  statistics,  the  transition 
function  is  time  homogeneous  and  we  have  for  a  grain-grain  transition 

P(At,  z,  E)  =  P{C :  E  3  £(At)  =  Cz)  (1.18) 

C  in  the  above  is  a  transmission  coefficient.  The  Huygens’  minor  principle  for  this 
process  is  expressed  by  a  semi-group  property  . 

W(t1+t2)  =  W(t1)W(t2)  (1.19) 


of  the  operator 


W(t)[f(z)]  =  J  f(z)P(t,z,  dz),  f  e  L!  (Z)  (1.20) 

z 


Also  here  we  introduce  a  diffusion  approximation.  This  is  expressed  by 
i)  the  Fokker-Planck  equation 

■^p(t,  z)  =  -  -|-[A;(z)p(t,  z)]  +  Y^-[B^(z)p(t,  z)]  (1.21) 


in  which 
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ii)  p(t,  2)  is  the  probability  density 


p(t,  z)  =  Jp(t,  z\  z)dz  (1.22) 

Z 

iii)  satisfying  the  initial  condition 

p(0,  z)  =  8(z  -  Zq)  (1.23) 

iv)  and  p(t,  z  ,  z)  being  such  that 

A;(z)  =  _!_<(<  “t>)  -  C(0)  |  C(0)  =  z>  (1.24) 

<“x> 

BW(z)  =  -J— <rt«“x»  -  m2 1  C(0)  =  z>  (1.25) 

<°T> 

-i—E<[C(<  °x>)  -  C(0)]n  |  C(0)  =  z>  =  0,  n  >  2  (1.26) 

<°t> 


Similar  to  (1.17),  also  (1.26)  represents  a  constraint  on  the  strength  of  randomness  of  the 
medium’s  microstructure.  Our  preliminary  calculations  indicate  that  (1.17)  is  more 
restrictive  than  (1.26).  These  calculations  were  based  on  the  assumption  that  the  random 
variables  E  and  p  have  independent  uniform  probability  densities. 

We  observe  from  (1.1)  that  in  case  of  a  nonlinear  constitutive  law,  the  2;t  process  is 
dependent  on  the  £t  process.  Thus,  is  then  Markov  only  in  the  sense  of  the  vector 
process  (£,  £)t  being  Markov.  In  this  case  we  have  the  transition  probability  function 

P(At,  (z,  x),  E,  F)  =  P(C(At)  €  E,  £(At)  e  F|C  =  z,$  =  x)  (1.27) 
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2.2  Wavefront  Propagation  Due  to  an  Arbitrary  Pulse:  Linear-Elastic  Grains 

For  the  case  of  a  microstructure  with  grains  whose  constitutive  law  is  piecewise 
linear,  the  process  is  not  explicitly  dependent  on  the  £t  process,  providing  £  takes 
values  in  any  given  linear  range;  we  consider  linear-elastic  grains  here.  Thus  both 
processes  are  Markovian  separately,  and  their  transition  functions  are  given  by  (1.4),  and 
(1.15),  respectively.  However,  their  joint  probability  density  p(t,  x,  z)  depends  on  the 
joint  transition  function  p(t,  x  ,  z,  x,  z)  as  expressed  by  the  conventional  relation 

p(t,  x,  z)  =  J  J  p(0,  x',  z')p(t,  x\  z',  x,  z)dx'dz'  (2.1) 

HZ 


where 


„  _\  _  d2P(ti  x  .  z\  x,  z) 
p(t,  x  ,  z ,  x,  z)  =  — - — - 


(2.2) 


and  p(0,  x  ,  z  )  is  the  initial  probability  density. 

We  now  consider  the  diffusion  approximation  for  p  =  p(t,  x,z).  The  Fokker-Planck 
equation  reads 


•|p=-|r[A5P]-|-[A5P] 


dx 


dz 


+  i^-tBuPi  +  +  y3sr[B«pl + 


(2.3) 


Clearly,  a  complete  diffusion  process  description  of  the  pulse  evolution  is  obtained  with 
(2.3)  rather  than  with  a  system  of  two  equations  (1.9)  and  (1.18)  -  this  is  due  to  the 
presence  of  "cross-terms,"  i.e.  fourth  and  fifth  on  the  right  hand  side  above.  In  the 
following  we  adopt,  for  simplicity,  the  approach  based  on  the  system  of  (1.9)  and  (1.18). 
As  a  matter  of  fact,  this  corresponds  to  following  the  evolution  of  the  pulse  in  terms  of  an 
operator 
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T(t)  =  U(t)W(t) 


(2.4) 


Considering  first  the  process  we  need  A^(x)  and  B^(x)  in  terms  of  the 
medium’s  statistics.  From  (1.8),  (1.15)  and  (1.16)  we  obtain* 

A5(x)  =  0,  Bs(x)  =  =  < A/f  »<> =  D  <2-5) 

where,  for  simplicity  of  notation,  we  have  introduced  a  constant  D.  Now,  (1.12)  becomes 
a  simple  diffusion  equation 

l-pft.x)  =  I  BS  p(t,x)  (2.6) 

Solution  to  (2.6)  subject  to  the  initial  condition 

p(t,x)  =  8(x)  at  t=0  (2.7) 


is,  as  is  well  known, 

l 

p(t,x)  =  (2jtB^t)  2  exp  [-x2/2B^t]  (2.8) 

Turning  now  to  the  process,  we  need  to  determine  the  coefficient  functions 
A^(z)  and  B^(z).  These  will  depend  on  the  transmission  coefficient  “C  and  the  relative 
impedance  11  £ 


*  «•»  denotes  variance 
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(2.9) 


UC  = 


2hX 

1+V 


In  the  above,  i  and  t  denote  the  regions  (i.e.  grains)  of  incident  and  transmitted  waves, 
respectively.  From  (1.18),  (1.24)  and  (1.25)  we  find 


A?(z)=  <i|>  (<Uc>“1)=Az 


(2.10) 


(<itC2>-2<i,C>+  l)  =  Bz2 


where,  for  simplicity  of  notation,  we  have  introduced  two  constants  A  and  B.  It  follows 
that  the  Fokker-Planck  equation  (1.21)  is,  in  Van  Kampen’s  (1981)  terminology,  a  so- 
called  nonlinear  equation;  the  cause  of  the  nonlinearity  lies  in  the  multiplicative  character 
of  the  process  as  reflected  by  the  transition  function  (1.18).  At  this  stage  we  can  either 
seek  the  solution  p(t,z),  or  assess  the  time  dependence  of  the  first  and  second  moments 
only.*  We  decide  here  on  the  latter  alternative,  and  hence  obtain  from  (1.21) 


=  <A^(z)> 

=  2<(C  -  <>)A;(z)>  +  <B%(z)> 


(2.11) 


These  have  to  be  supplemented  with  the  initial  conditions 


*  Solution  of  the  equation  (2.3)  is  obtained  in  Ostoja-Starzewski  (1991). 
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<£>  =  Zq  and  <:£:*>  =  ()  at  t=0 


(2.12) 


In  view  of  (2.10),  relations  (2.11)  become 

=  (2A  +  B)<C»  +  B  <>2 

Finally,  we  obtain  (note  from  (2.12)  that  <C(0)»  =  0) 

<C(0>  =  Zo  exp  (At) 


(2.13) 


(2.14) 


In  this  section  we  study  wavefront  propagation  in  a  microstructure  made  of  linear- 
elastic  grains,  occupying  the  X  k  0  half-space,  and  subjected  to  a  surface  pressure  f(t)  at 
X  =  0,  see  Fig.  2.  The  constitutive  law  of  any  grain  is 

<x  =  E(co)e  (2.15) 

where  (0  indicates  that  the  Young’s  modulus  E  is  a  random  variable.  The  above 
statement  corresponds  to  a  problem  of  wave  propagation  in  a  1-D  rod  made  of  elastic 
elements  with  varying  moduli.  By  introducing  the  constitutive  law 

o  =  [A<to)  +  2n(co)]e  (2. 16) 

we  have  a  problem  of  wave  propagation  in  a  half-space  made  of  layers  of  varying 
(random)  moduli  X  +  2\l  For  simplicity  of  notation  we  stay  with  equation  (2. 15). 

There  are  two  regions,  I  and  II,  in  the  space-time  domain  (Fig.  3).  In  region  I  we 
have  for  stress,  strain,  and  particle  velocity: 
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Oi  =  0,  £i  =  0,  Vj  =0 


(2.17) 


For  any  specific  medium  B(o)),  region  I  is  bounded  by  a  trajectory  above  which  the  body 
is  disturbed;  the  latter  region  is  denoted  by  n.  The  stress,  strain,  and  particle  velocity  in 
region  II  are  determined  jointly  by  the  ^-process,  the  constitutive  law  (2.15)  and  the 
well-known  relation 


v  =  —  c£  (2.18) 

If  the  pressure  is  prescribed  it  is  most  natural  to  represent  stress  by  £.  In  view  of  the 
developments  of  section  2.1,  we  obtain  a  family  of  boundaries  X(0,  co)  separating  both 
regions,  that  is,  undisturbed  I  from  disturbed  IL  Rather  than  working  with  this  "fuzzy" 
ever-broadening  boundary  we  introduce  the  mean  path  OA  and  trace  the  character  of 
evolution  of  the  disturbance  in  the  random  medium  B  =  {B(to);  co  e  ft}  in  terms  of  the 
process  referred  to  this  path. 

Now,  it  follows  by  the  argument  of  linearity  of  response  of  all  the  grains  (see 
equation  (2.15))  that  the  propagation  speed  c(0>)  is  same  for  every  stress  level  in  any 
given  grain.  Hence,  for  any  fixed  time  t,  the  trajectory  X(t,  to)  is  same  as  X(0,  co)  for  any 
given  body  B(co),  except  for  the  time  shift  t.  This  is  analogous  to  two  characteristics 
being  parallel  in  case  of  wave  propagation  in  a  homogeneous  deterministic 
microstructure.  In  fact,  it  follows  that  the  reference  paths  are  parallel  for  a  linear  elastic 
microstructure. 

2.3  Wavefront  Propagation  Due  to  a  Square  Pulse:  Linear-Hysteretic  Grains 

In  this  section  we  study  wave  propagation  in  a  body  whose  elements  (i.e.  grains  in 
case  of  a  rod  model,  or  layers  in  case  of  a  half-space  model)  are  governed  by  a  stress- 
strain  diagram  of  Fig.  4.  The  stress-strain  curve  is  a  straight  line  OM  on  initial  loading  to 
M;  its  slope  defines  the  initial  modulus  E0.  Upon  unloading  the  stress-strain  curve  is 
another  straight  line  MN  which  defines  the  second  modulus  Hj 


15 


O  —  Om  =  Ej  (£  —  £,„) 


(3.1) 


If  material  is  reloaded,  it  follows  the  line  NM  to  M,  and  then  continues  along  the  initial 
loading  line.  E0  and  Ej  are  random  from  grain  to  grain;  they  correspond  to  phase 
velocities  Co  and  ci,  respectively. 

We  consider  wave  propagation  under  action  of  a  pressure  pulse  f(t)  applied  at 
X  =  0,  where  f(t)  is  a  square  pulse 


f(t)  = 


Po 

0 


for  0  £  t  <  ti 
otherwise 


(3.2) 


The  problem  with  this  initial  condition  for  a  deterministic  homogeneous  medium,  which 
was  solved  by  Salvador!  et  al  (1960)  (see  also  Rohani,  1970),  forms  the  reference  basis 
for  solution  of  a  stochastic  problem.  In  Fig.  5  we  give  the  space-time  diagram  of  the 
deterministic  problem.  We  see  that  there  are  several  regions:  I,  n,  m,  IV,  V, ....  Region 
I  is  that  of  an  undisturbed  body,  i.e. 


01=0,  ei  =0,  V!  =0  (3.3) 

Region  n  corresponds  to  a  material  in  which 

°n  =  ”  Po*  en  =  oh/Eq,  vjj  = — c0en  (3.4) 

the  boundary  between  I  and  II  being  given  by 

X  =  Cot  (3.5) 


In  region  HI  we  have 


On  =  0,  vra  =  Vj  - 


Po 

pcj 


(3.6) 


The  time  t2  at  which  the  pressure  p,  reaches  the  original  wavefront  is  obtained  as 
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Cl 


(3.7) 


Cl  -C0 


tl 


As  discussed  in  the  above  reference,  at  time  t2  there  is  a  reflection  of  the  wavefront 
travelling  at  velocity  ci  (i.e.  upper  boundary  of  region  II)  from  the  original  wavefront 
travelling  at  ce  so  that  a  new  region  IV  is  created.  The  field  quantities  in  IV  are 


Op/  =  apo,  vjv  =  -  oiv/pci 


(3.8) 


where 


Ci/c0  - 1 
ci/c0  + 1 


(3.9) 


The  analysis  can  be  continued  in  this  fashion  to  cover  the  entire  X,  t-plane,  and  to 
study  arbitrary  pulses  f(t). 

Turning  now  to  the  random  medium  problem  we  see  that  each  of  the  lines  in  Fig.  5 
representing  the  discontinuity  waves  (i.e.  shocks)  can  be  considered  as  a  mean  path 
providing  reference  for  stochastic  processes  £  and  £.  Thus,  for  example,  the  leading 
shock  is  a  reference  for  a  family  of  characteristics 

X(0,  ca,  +  c0)  =  X(t,  G> ,  4c0) | t=o>  coe  Q  (3.10) 

Similarly,  the  line  bounding  the  region  II  from  above  is  a  family  of  characteristics 

X(t1,a,+cl)  =  X(t,©,+c1)|t=tl,©e  Q  (3.11) 

and  so  forth.  The  cones  corresponding  to  the  %  processes  along  the  above  mentioned  two 
mean  characteristics  (paths)  are  shown  in  Fig.  6.  Clearly,  the  point  of  intersection  will  be 
diffused  about  the  reference  point  (X  =  <c0>t2,  t2). 

First,  we  consider  one  characteristic  X(0,  to,  +c0)  which  is  slower  than  the  mean 
one.  This,  we  may  say,  corresponds  to  an  effective  velocity  c0  =  c0tef  for  body  B(co) 
being  smaller  than  <c0>;  <c0>  defines  the  average  path  according  to  (1.2).  In  case  the 
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moduli  Eo(<0)  and  Ei(co)  in  the  constitutive  law  are  independent  random  variables,  there 
is  no  way  to  predict  whether  the  characteristic  X(tj ,  co,  +Cj )  for  the  same  body  B(co)  will 
also  be  slower  than  the  mean  one  corresponding  to  <+ci>.  In  that  case  the  intersection  of 
both  characteristics  -  X(0,  co,  +c0)  and  X(tlt  co,  +Cj)  -  may  take  place  anywhere  in  the 
quadrilateral  ABFE.  In  the  special  case  of  dependence  between  random  variables  E0(co) 
and  Ei  (co)  being  such  that 

Eo(co)  <  <Eo>  if  and  only  if  Ei(co)<<Ei>  (3.12) 

and  the  same  holding  for  ">H  signs  in  place  of  the  point  of  intersection  will  fall  in  the 
HBFO'  quadrilateral. 

Now,  if  we  consider  a  characteristic  X(0,  co,  +cD)  which  is  faster  than  the  mean 
one,  we  obtain  immediately  the  following  conclusions  in  a  complete  analogy  to  the 
above  analysis: 

i)  for  no  dependence  of  moduli  Eo(co)  and  Ei(co),  the  intersection  will  occur 
somewhere  in  the  quadrilateral  EDCF, 

ii)  for  a  dependence  of  moduli  according  to  (3.12),  the  intersection  will  occur 
somewhere  in  the  quadrilateral  EDGO  . 

It  is  seen  from  the  geometry  of  the  problem  that,  in  case  of  dependence  between 
Eo(co)  and  Ei(co)  for  all  co,  the  scatter  in  point  of  intersection  of  two  characteristics  -  as 
measured  by  the  distance  between  points  AC  -  increases  for  the  ratio  E}/E0  decreasing  to 
1.  This  indicates  that  even  a  weak  randomness  in  the  medium’s  properties  may  alter 
certain  aspects  of  its  response  in  a  significant  way. 

On  the  other  hand,  in  case  of  dependence  defined  by  (3.12),  the  scatter  in  the 
location  of  point  of  intersection  of  two  characteristics,  as  measured  by  the  distance 
between  points  B  and  D,  is  much  weaker. 

The  foregoing  analysis  of  the  intersection  of  characteristics  for  a  random 
microstructure  governed  by  a  bilinear-type  constitutive  law  carries  over  to  any 


18 


intersection  point  of  Fig.  5,  as  well  as  to  all  such  points  of  other  bilinear-type  (and 
piecewise-linear)  problems. 

2.4  Numerical  Results 

In  Section  2.2  we  have  derived  the  diffusion  equations  governing  the  behavior  of 
probability  densities  p(tfx)  and  p(t,z).  From  these  we  have  obtained  explicit  formulas  for 
p(t,x)  and  the  first  and  second  moments  of  the  variable  £,  that  is  equations  (2.8)  and 
(2.14).  These  formulas  are  governed  by  the  coefficients  D,  A,  and  B,  which,  in  turn, 
depend  on  the  three  random  variables  1,  p,  and  E  describing  the  microstructure.  It  is  of 
interest  now  to  investigate  the  effects  of  this  dependence. 

Consonant  with  our  assumption  of  a  strict  spatial  homogeneity  of  material 
properties  (recall  the  discussion  preceding  equation  (1.4))  we  describe  the  properties  of 
any  grain  by  a  density  p(l,p,E),  and,  in  what  follows,  assume  a  particular  special  case 

p(l,p,E)  =  p(I)p(p)p(E)  (4.1) 

Calculations  of  D,  A  and  B  according  to  (2.5)2  and  (2.10)  are  unwieldy.  Indeed, 
even  the  calculation  of  the  probability  density  of  the  phase  velocity  c  =  ^E/p  is  lengthy. 
The  problem  qualifies  perfectly  for  an  application  of  a  Monte  Carlo  technique.  Results 
based  on  this  approach  are  summarized  in  Tables  1  and  2;  C  denotes  <c>  here.  These 
calculations  correspond  to  the  situation  of  p(l),  p(p),  and  p(E)  being  uniform  densities 
with  means  equal  1  and  standard  deviations  Oi,  ap,  and  Oe  ranging  from  0.02  to  0.2.  In 
order  to  bring  out  the  relative  importance  of  randomness  of  1  versus  p  versus  E,  four 
general  cases  were  considered: 

i)  d  =  cp  =  Oe  =  0.02  through  0.2:  left  part  of  Table  1 , 

ii)  G\  -  0,  op  =  Oe  =  0.02  through  0.2:  right  part  of  Table  1 , 
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iii)  Oj  =  Og  =  0,  ap  =  0.02  through  0.2:  left  part  of  Table  2, 


iv)  ci  =  ap  =  0,  Oe  =  0.02  through  0.2:  right  part  of  Table  2. 


TABLE  1 


0[  — Op  — Oe 

A 

B 

c 

D 

<Ji=0 

<*p=<7E 

A 

B 

C 

D 

.0005 

1.0002 

.0002 

0.02 

.0005 

.0002 

.0002 

1.0003 

.0008 

0.04 

.0002 

-.0009 

.0005 

1.0009 

.0019 

0.06 

.0005 

1.0010 

.0018 

0.08 

-.0014 

.0008 

1.0016 

.0033 

0.08 

-.0014 

.0013 

1.0024 

.0050 

0.10 

-.0019 

.0012 

1 

■  7= 

0.12 

-.0044 

.0019 

1.0043 

.0075 

0.12 

.0018 

■mm 

-.0053 

1.0049 

.0102 

0.14 

.0025 

.0100 

WSM 

-.0075 

7 

1.0095 

.0135 

0.16 

.0032 

.0134 

mm 

-.0097 

1.0076 

.0171 

0.18 

.0043 

0.20 

1.0092 

0.20 

-.0102 

.0052 

1.0132 

.0219 

TABLE  2 


d=cE=0 

op 

A 

B 

C 

D 

<*E 

A 

B 

c 

-.0000 

I 

.0001 

-.0000 

.0000 

1.0000 

.0000 

-.0002 

.0004 

-.0002 

.0001 

.9998 

.0004 

0.06 

-.0004 

1.0013 

.0009 

0.06 

-.0004 

.0002 

.9996 

.0009 

0.08 

-.0011 

.0004 

1.0030 

.0016 

0.08 

-.0011 

.0004 

.9986 

.0016 

0.10 

-.0012 

.0006 

1.0038 

.0026 

0.10 

-.0013 

.0006 

.9988 

.0025 

0.12 

-.0018 

.0009 

1.0054 

.0037 

0.12 

-.0018 

.0009 

.9983 

.0037 

0.14 

-.0024 

.0013 

1.0077 

.0052 

0.14 

-.0025 

.0013 

.9978 

.0050 

0.16 

-.0028 

.0017 

1.0089 

.0068 

0.16 

-.0027 

.0017 

.9974 

.0065 

0.18 

-.0036 

.0021 

1.0117 

.0089 

0.18 

-.0036 

.0021 

.9969 

.0083 

0.20 

-.0047 

.0027 

1.0149 

.0113 

0.20 

-.0046 

.0027 

.9959 

.0104 

The  numbers  of  both  tables  were  obtained  by  using  104  samples  for  each  specific 
case  (i.e.  each  specification  of  all  three  c’s).  The  last  digits  are  only  approximate  since  a 
rather  simple  random  number  generator  was  employed. 

The  results  of  Table  1  show  that  the  effect  of  random  grain  lengths  vis-k-vis  the 
deterministic  ones  is  practically  negligible.  On  the  other  hand,  restricting  the  material 
randomness  to  p  alone,  or  E  alone,  has  a  definite  effect  of  decreasing  the  randomness  of 
the  system:  weaker  attenuation  and  smaller  diffusion.  Interestingly,  case  iv)  is  the  only 
one  which  results  in  <c>  <  1. 

The  above  being  a  dimensionless  formulation  it  clearly  offers  a  possibility  of 
calculating  any  particular  dimensional  case  if  the  means  p  and  standard  deviations  o  of  1, 
p,  and  E  are  given.  Thus,  for  example,  if  the  following  case  is  given: 

p,  =  l"  =  0.0254m  ,  Gi  =  0.0 

Up  =  100  lb/ft3  =  1602  ap  =  0.1  •  Pp  (4.2) 

m 

He  —  63  ksi  =  3016440— L-,  Oe  —  0.1  •  Pe 

nr 

we  read  from  the  right  part  of  Table  1  the  following  values: 

A  =  -0.0019,  B  =  0.0012,  C  =  1.0023,  D  =  0.0052  (4.3) 

V  V  V  V 

Now,  the  dimensional  values  of  the  four  constants  -  denoted  by  A,  B,  C,  and  D  -  can  be 
obtained  from  these  formulas 
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(4.4) 


Thus,  (4.3)  yields  (in  the  SI  system) 


V 


A  = 


v  1 
B  =  2.05— 
s 


(4.5) 


C  =  43.49— 
s 


D  =  9.79 


nr 
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m.  FRAMEWORK  FOR  THE  WAVEFRONT  PROPAGATION  STUDIES  IN  2-D 
AND  3-D  RANDOM  MICROSTRUCTURES 

3.1  Basic  Model 

The  stochastic  model  of  disturbance  evolution  presented  in  the  preceding  section 
forms  a  basis  for  a  treatment  of  wavefronts  propagating  in  2-D  and  3-D  media.  For 
illustration  of  the  analysis,  we  consider  the  2-D  case  here  only.  Specifically,  we  look  at 
microstructures  representable  by  graphs,  such  as  the  one  shown  in  Fig.  la).  As  discussed 
by  Ostoja-Starzewski  (1987,  1989a),  thin  lines  connect  the  centers  (vertices  of  the  set  V) 
of  interacting  grains  and  form  set  E  of  graph  G(V,  E).  On  the  other  hand,  thick  lines  of 
the  set  E  outline  the  contours  of  polygon-shaped  grains;  they  meet  in  triplets  at  vertices 
of  the  set  V ,  and  hence  we  have  the  graph  G(V  ,  E ).  It  is  assumed  that  the  space  of  any 
polygon  domain  is  filled  by  a  homogeneous  continuum  so  that  there  are  no  voids.*  The 
continuum  itself  is,  in  general,  anisotropic  and  hence  the  propagation  characteristics  of 
any  type  of  wave  (e.g.,  dilatational)  are  given  by  an  envelope  of  velocity  vectors,  i.e.  an 
indicatrix,  specific  to  any  grain.  Thus,  the  microstructure  can  be  characterized  by  a  set  of 
indicatrices,  each  of  which  is  attached  at  the  vertex  of  G(V,  E),  see  Fig.  lb). 

The  microstructure  represented  in  Fig.  la)  and  lb)  is  just  one  possible  realization 
of  microscale  geometrical  and  physical  properties  of  the  medium  in  this  part  of  spatial 
domain,  and  therefore  it  is  a  deterministic  medium  B(co);  co  indicates  a  parameterization 
by  an  element  from  the  sample  space  £2.  Accordingly,  a  random  medium  B  is  a  family 
(B(co);  toe  Q)  of  all  bodies  thus  described,  with  probability  distributions  of 
microstructural  geometries  and  physical  characteristics  specified.  In  the  following  we 

*  This  model  permits  a  simulation  of  actual  voids  (i.e.  absence  of  grains)  by  letting  the  elastic  moduli  of  a 
given  polygon  domain  be  equal  to  zero. 
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assume  these  probability  distributions  to  be  spatially  homogeneous. 

Evolution  of  disturbances  follows  the  kinematics  of  rays  governed  by  two 
equations 


dX  =  c2Ydt 


dY  =cV 


(1.1) 


where  X  and  Y  are  the  position  and  direction,  respectively,  of  the  ray.  Evidently,  the 
pair  (X,  Y)  undergoes  a  Markovian  evolution,  and  we  can  write  its  transition  function  in 
terms  of  the  medium’s  statistics.  However,  it  is  more  convenient  to  formulate  this 
process  in  terms  of  the  dispersion  vector  (recall  (1.2)) 

£(t,co)  =  X(t,  co)-<X(t)>  (1.2) 

and  the  angle  tj  shown  in  Fig.  7.  In  this  figure  we  give  a  space-time  graph  of  plane 
disturbance  propagation  in  Xi,  X2 -space,  where  the  initial  conditions  are 

X(0)  =  (0,  0),  Y(0)  =  1(1,0). 

~  ~  c 

Introducing  the  transition  probability  function 

P(t,x,y,E)  =  P{[|(t),Tl(t)]e  E||(0)  =  x,Tl(0)  =  y}  (1.3) 

we  have  a  Markovian  operator 

U(t)[f(x\  y')]  =  J  J  f(x,  y)P(t,  x',  y',  dx',  dy'),  f  e  Lx  (E  x  H)  (1.4) 

_H  - 


which  satisfies  the  semi-group  property 


U(t1+t2)  =  U(t1)U(t2) 


(1.5) 
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Similar  to  the  1-D  model,  the  spatial  inhomogeneity  of  the  medium  has,  besides 
randomization  of  characteristics,  an  effect  on  wave  amplitude  carried  by  any  disturbance. 
In  Fig.  8a)  we  show  a  typical  grain  a  of  the  microstructure  in  which  a  plane  disturbance 
propagates.  The  modulation  of  the  disturbance  amplitude  in  its  passage  through  grain  a 
is  described  by  a  transmission  operator 

C:  “£(0-*  “C(t  +  At)  (1.6) 

Here,  again,  £  is  a  field  quantity  such  as  stress  or  strain,  and  i  and  t  indicate  the  incident 
and  transmitted  quantities,  respectively.  Evolution  of  the  disturbance  takes  place  along 
the  actual  random  ray  path  specific  to  the  grain  a;  the  average  path  is  shown  for 
reference.  The  development  at  the  end  of  Section  2.1  applies  here  as  well,  and  thus  W(t) 
is  the  appropriate  evolution  operator. 

3.2  Spherical  and  Cylindrical  Waves,  and  Local  Averaging 

The  approach  presented  in  the  1-D  setting  (Section  II)  may  be  briefly  described  as 
a  combination  of  the  method  of  characteristics  with  the  diffusion  process  approximation 
of  arrival  times  and  field  quantities  referred  to  the  mean  paths  (i.e.  average 
characteristics).  This  suggests  a  method  for  2-D  and  3-D  settings:  use  the  solutions  to 
deterministic  wave  problems  with  spherical  or  cylindrical  symmetry  as  a  reference  for 
the  diffusion  processes  (£,T])t  and  £t  in  order  to  assess  the  results  (e.g.  strength  of 
fluctuations)  in  random  media  (books  by  Nowacki  (1978)  and  Wlodarczyk  (1986)  are 
good  references  on  such  problems).  In  these  problems,  on  account  of  the  assumed 
symmetry,  all  the  field  quantities  are  functions  of  time  and  a  single  spatial  variable 
(radius  r).  However,  higher  level  of  complexity  is  encountered  here  due  to  the  presence 
of  more  complex  stress  states  than  those  in  the  1-D  problems.  Of  course,  any  type  of 
spatial  symmetry  is  not  a  condition  for  the  applicability  of  our  Markovian  theory,  but 
then  more  complicated  analyses  need  to  be  carried  out. 


25 


Finally,  we  note  that  the  entire  approach  adopted  here  gives  descriptions  of  the 
wavefront  (t)  -  random  field  in  space-time  -  with  the  scale  of  resolution  equal  to  the 
average  grain  size  <1>.  Evolution  of  4^(0  is  described  along  any  given  average 
characteristic  by  the  propagator  T(t)  such  that  (Huygens’  minor  principle) 

T(ti+t2)  =  T(ti)T(t2)  (1.7) 

This  corresponds,  in  fact,  to  a  local  averaging  of  the  true  wave  process  at  the  scale  <1>. 
By  using  the  notion  of  local  averaging  of  random  fields  (Vanmarcke,  1983)  we  can  arrive 
at  less  detailed  but  smoother  descriptions  of  0(0,  all  of  them  parameterized  by  8  =  L/<1>, 
where  L  is  the  scale  of  averaging.  The  situation  is  akin  to  the  notion  of  a  window  in  die 
problem  of  continuous  parameter  random  field  approximations  in  quasi-statics  of  random 
microstructures;  see  Ostoja-Starzewski  and  Wang  (1989a  and  b).  In  the  latter  case  L  is 
the  edge  length  of  a  window  AV,  AV  being  a  representative  volume  element.  Clearly,  for 
any  8  <  <»  we  have  a  statistical  continuum,  while  for  5  — ►  we  obtain  a  deterministic 
continuum  since  the  number  of  microelements  (e.g.  grains)  in  AV  becomes  infinite. 
Returning  to  our  present  problem  we  show  a  ray  tube  of  width  L  or  5  in  the 
microstructure  of  Fig.  8b).  Hence,  the  analysis  can  be  conducted  in  terms  of  a 
transmission  operator  C(8)  and  a  propagator  T(t,  8)  corresponding  to  a  given 
approximation  8. 


26 


IV.  CONCLUSIONS 


In  this  report  we  outlined  a  model  for  analysis  of  stress  wavefronts  propagating  in 
granular  media  without  voids.  The  following  are  two  keystones  of  the  model: 

—  microstructure  is  representable  by  a  graph, 

—  attention  is  focused  on  the  evolution  of  a  single  point  of  wavefront  (pulse)  called  a 
disturbance. 

The  report  is  based  on  our  previous  work  on  wave  propagation  in  discrete  random 
media.  Specifically,  we  use  the  result  that  the  evolution  of  a  forward  propagating 
disturbance,  travelling  along  its  stochastic  characteristic,  is  described  by  a  random  vector 
process  (£,£) t;  £  is  the  amplitude  of  the  field  quantity  (e.g.  stress)  at  the  disturbance,  and 
^  is  the  fluctuation  in  distance  covered  by  the  disturbance  up  to  time  L  Thus,  with  £  we 
model  stochastic  temporal  variability  of  the  amplitude  of,  say,  stress,  while  with  £  we 
model  diffusion  of  characteristics  in  space-time,  and  hence,  scatter  in  the  arrival  times. 

(£,£)t  is  a  Markov  process  for  a  microstructurc  with  grains  of  a  nonlinear  type, 
where  nonlinearity  reflects  either  a  non-Hookean  nature  of  elastic  grains,  or  a  nonelastic 
nature  (e.g.  elastic-plastic)  of  their  constitutive  laws.  In  the  case  of  a  nonlinear 
constitutive  law  there  arises  an  interesting  phenomenon  of  curving  of  the  effective 
average  paths  of  forward  propagating  waves;  this  is  presently  under  study  and  will  be 
reported  in  the  next  paper.  It  is  important  to  note  that,  (£,£),  is,  strictly  speaking,  Markov 
only  for  grains  whose  properties  are  not  mutually  dependent  in  space,  i.e.  when  the 
physical  properties  of  the  grains  (mass  density  and  constitutive  moduli)  are  independent 
random  variables.  Otherwise,  the  process  (£.5)i  retains  its  Markov  property  only  in  the 
sense  of  being  driven  by  the  random  Markov  process  of  physical  properties. 

A  very  practical  approximation  of  the  (£,£)t  process  is  obtained  for  grains  whose 
constitutive  law  is  piecewise  linear.  The  simplest  approximation  is  obtained  in  the  form 
of  two  Fokker-Planck  equations  -  one  for  and  one  for  -  in  which  the  key  coefficient 
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functions  can  be  calculated  explicitly  in  terms  of  the  statistics  of  medium’s  properties. 
The  process  is  then  governed  by  a  linear  diffusion  equation,  and  hence  its  probability 
density  at  any  arbitrary  time  t  >  0  is  Gaussian.  The  £t  process  is  a  multiplicative  one  -  its 
solution  is  found  here  in  terms  of  the  first  and  second  moments  only;  complete  solution 
for  the  probability  density  is  possible  too.  The  rates  of  change  of  field  quantities  in  the 
model  are  referred  to  the  temporal  microscale  in  the  problem:  At  =  <l/c>. 

After  discussing  the  implications  of  the  above  stochastic  model  for  a 
microstructure  of  linear  elastic  grains,  we  turn  to  the  case  of  linear-hysteretic  grains.  We 
develop  a  method  of  solution,  also  applicable  to  all  transient  wave  problems  with 
microstructurcs  of  piecewise  linear  constitutive  laws,  which  uses  the  space-tune  diagram 
of  the  deterministic  medium  problem  as  a  reference  for  the  average  paths  of  the 
stochastic  setting.  It  is  shown  that  there  can  occur  a  very  strong  scatter  in  solutions  due 
to  an  even  small  randomness  of  the  medium. 

In  the  presentation  of  our  recent  conference  paper  (Ostoja-Starzewsld,  1989b)  we 
have  discussed  an  extension  of  this  analysis  to  microstructures  governed  by  bilinear- 
elastic  and  nonlinear  elastic  laws  with  random  constitutive  coefficients.  In  the  first  case 
there  is  a  possibility  of  an  abrupt  change  in  the  orientation  of  the  forward  characteristics, 
while  in  the  second  case  there  is  a  continuous  curving  of  characteristics;  see  (Ostoja- 
Starzewsld,  1991)  for  analysis  of  interesting  phenomena  which  arise  in  the  first  situation. 

In  a  separate  section  of  the  report  we  present  results  of  calculation  of  the  key 
coefficients  (drift  and  diffusion)  which  appear  in  the  equations  governing  the  evolution  of 
forward  characteristics  and  modulation  of  pulse  strength.  The  calculations  are  carried 
out  in  the  dimensionless  setting  of  a  1-D  linear  elastic  model.  This  has  an  obvious 
advantage  that  any  physical  (dimensional)  case  can  be  calculated  with  the  formulas  (4.4) 
and  Tables  1  and  2,  under  the  condition  that  the  actual  case  satisfies  the  assumptions  of 
statistical  independence  of  generic  random  variables  1,  p,  and  E  describing  the 
microstructure. 
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In  the  final  part  of  the  report  we  outline  a  framework  for  solving  2-D  and  3-D 
problems.  First,  we  present  a  generalization  of  the  disturbance  propagation  process  in 
terms  of  vector  Markov  processes.  Analogous  to  the  1-D  setting,  this  can  be  combined 
with  the  classical  solutions  of  deterministic  problems  -  typically  well  known  in  cases  of 
spherical  and  cylindrical  symmetry  -  to  obtain  solutions  of  stochastic  problems.  Finally, 
we  briefly  introduce  the  concept  of  local  averaging  of  wavefronts  ~  an  idea  which  should 
be  very  useful  in  matching  experimental  measurements  with  this  theory. 

It  is  our  opinion  that  a  model  of  transient  wave  propagation  in  randomly 
heterogeneous  materials  has  to  be  derived  from  micromechanical  considerations.  The 
key  observation  is  that  a  propagating  disturbance  is  effected  by  the  random 
microstructure  of  the  material.  Hence,  one  has  to  be  very  careful  with  applying  effective 
continuum  models  that  were  initially  derived  for  quasistatic  situations.  As  the  analysis  in 
this  report  shows  the  random  nature  of  micro-scale  properties,  even  with  a  weak  strength 
of  fluctuations,  may  result  in  rather  strong  macroscopic  effects. 
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a.  Graph  representation  of  the  microstructure 
of  a  heterogeneous  solid 


Figure  1 


b.  Distribution  of  indicatrix  envelopes 
for  the  same  ensemble  of  grains 


Figure  2.  Space-time  graph  of  1-D  distribution  propagation 


Figure  5.  Space-time  graph  of  response  of  a  linear-hysteretic 
medium  to  a  square  pulse 


Figure  6.  Intersection  of  forward  dependence  cones  in  response  of 
a  random  linear-hysteretic  medium  to  a  square  pulse 
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Figure  7.  Space-time  graph  of  2-D  distribution  propagation 


38 


a.  Wave  passage  through  a  grain 


b.  Wave  passage  through  a  microstructure 
Figure  8 
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